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PERTURBATION OF DOMAIN: SINGULAR RIEMANNIAN 

METRICS 

C. MASON 

o. 

f^ , Abstract. We develop further some aspects of the spectral theory of a class 

OQ ' of Riemannian manifolds introduced by E. B. Davies; in particular we study the 

best constant in the Hardy Inequality which has become important in spectral 
theory. One application of this constant has been to a certain type of domain 
perturbation. This technique is useful when the domain has an irregular bound- 
ary and this is the case with the manifolds under consideration. However, in this 
paper we show that the manifolds possess a Hardy constant that lies outside the 
Ph ' range permitted by existing theorems. Yet we are still able to prove theorems 

^^ I which give information about the domain perturbation problem and moreover, 

f~| ' we set up a specific example which can be used to show that our results are the 

jrt ! best possible. 



(N, 

^'- 1. Introduction 

00 

In the paper Two- Dimensional Riemannian Manifolds with Fractal Boundary, [|1|, 

^ . Davies introduced a class of Riemannian manifolds which exhibited unusual prop- 

(^ ! erties. These manifolds are simple in that the sense that there is only one coor- 

O I dinate patch and the geometry when viewed with the natural Euclidean metric 

^ ■ is extremely straightforward (and much of the corresponding spectral theory well 

c^ . understood). However, when equipped with a certain Riemannian metric, which 

becomes singular at the boundary, the manifolds show a link with Euclidean do- 
mains that have fractal boundaries. More precisely Davies constructs a piecewise 
smooth homeomorphism between the unit disc with a specific example of the sin- 
gular metric and the von Koch snowfiake domain equipped with the Euclidean 



X 

Cu metric. 



One strategy for computing the Dirichlet eigenvalues of a region with irregular 
boundary is to use the following perturbation of domain argument. Given a Rie- 
mannian manifold Q one considers, for small enough e > 0, a smaller domain Q^ 
such that 

{x e fi : dist(a;, dVl) > e} cn^cn. 

We denote by H^ and i^e.z? the Dirichlet Laplacians associated with Q and Q^ 
respectively. If Xn,e denotes the nth eigenvalue of H^^d and A^ denotes the nth 
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2 C. MASON 

eigenvalue of Hq then variational arguments imply that 

for all n. Our goal is to find estimates of the quantity Xn,e — A„ to assess how 
quickly these approximations converge as e -^ 0+. 



Results relating to this problem are given by Davies in and by Pang in |T2| 
and |TB|. Pang's results are applicable to the bottom eigenvalue for open, bounded 
and simply connected sets. Davies' results are valid for all eigenvalues and it is 
the techniques in this paper that we will follow (Pang's proofs are completely 
unrelated). The key assumption made in [^ is that a Hardy inequality of the type 

holds in the sense of quadratic forms. Here d{x) :=dist(x, 9f2) and c,a are some 
constants. The infinimum of all possible values for c such that there exists an 
a < oo for which the above inequality holds is known as the Hardy constant. 

The key result in is the following: 

2 

(1) A„,e - A„ < CnS' 

where c„ is some constant and c is the Hardy constant associated with the region 
which is assumed to satisfy 

c> 2. 



The restriction on c is important for the proof presented in [Q]. In Theorem |4.5| we 
prove that the class of singular manifolds introduced in ||I| has a Hardy constant 
that satisfies 

1 < c< 2. 

This result stands in contrast to the situation of regular Euclidean domains. It is 
known that the Hardy constant for every Euclidean domain whose boundary has 
at least one regular point is never smaller than 2. A precise formulation of this 
statement and proof may be found in Q] . 

Despite this we prove in Theorem ^^ that for any manifold with Hardy constant 
1 < c inequality ([|) is still true. The key difference is the use of better quadratic 
form and operator norm inequalities in Theorem |6]^. The conclusions as applied 
to the particular Riemannian manifolds are stated explicitly in Corollary |6.10| . We 
emphasise that the rate of convergence is faster than 0{e) when c < 2. 

Furthermore, in section 7 we give a concrete example. This is a rotationally in- 
variant domain and after separating variables one can reduce to a one dimensional 
problem. The techniques for dealing with the class of problems that arises are 
presented in another paper, [Q, since they are based on techniques from ordinary 
differential equations and are very different from anything presented here. They 
can be used to show that the constant achieved in our general method is sharp for 
this particular example. 
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The proofs require two distinct notions of the distance to the boundary and the 
relationship between them. The first distance is the usual Euclidean distance to 
the boundary and this is denoted by a. Since the boundary will be smooth in the 
Euclidean sense this is simply the length of the shortest straight line that intersects 
the boundary orthogonally. This function is only guaranteed to be smooth in a 
neighbourhood of the boundary but in general it satisfies 

in the distributional sense (we use the subscript here to stress that everything is 
the usual Euchdean quantity). 

The second distance is the distance induced by the Riemannian metric. This 
will be denoted d. This may be defined as the infimum of the lengths of paths 
connecting the point to the bounday where the length is calculated with respect 
to the Riemannian arc length. Theorem |3.3| gives a simple relationship between 
these two and is heavily used. 

2. The Basic Model 

We begin with Vt an open, bounded and connected subset of M.^ that has a C^ 
boundary and assign to it a Riemannian metric that becomes singular at the bound- 
ary in some controlled fashion. In future this will be referred to as the Riemannian 
metric. In terms of the arc length element and depending on some constant 7 to 
be specified later this is given by 

ds^ = a~''^^{x){dxl + ■ ■ ■ + dx%) 

where 

a{x) := distE{x,dfl). 

Remark 2.1. In the above and in all that follows a subscript E will denote a 
quantity taken with respect to the Euclidean metric. The absence of such a sub- 
script means that the Riemannian metric is appropriate. The fact that there are 
two metrics on the same domain is one of the principal challenges in the analysis. 

2.1. The Laplacian. The Riemannian volume element is given by 

dvol = a-^^{x)d^x. 
The constant 7 is chosen to satisfy the relation 

< A^7 < 1. 

Q therefore has finite volume and is incomplete. Moreover, the L^ spaces are 
defined to be the spaces for which the following norms are finite: 



p . 
p ■ 

'n Jn 



\f\Pdvol = / |/|V-^M^x. 



The norm on the tangent spaces derived from the Riemannian metric is such that 

|V/| = a^|Vs/U. 
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Since we have a metric we also have a Laplacian and this can be defined using 
quadratic form techniques. The quadratic form is given by 

Qif) := / Wfl'dvol = f WEflW-^'^^d'^x 
Jn Jn 

whenever this is finite. 

Definition 2.2. We define W^''^{Q) to be the space of weakly differentiable func- 
tions on Q for which 

Q(/) + ll/ll2<00 

and define W^'^n) to be the closure of C^{n) in W^'^Vt). 

Remark 2.3. If 7 = then these are the standard Sobolev spaces which we will 
denote by W^'^n) and W^^^i^) 

Definition 2.4. The Dirichlet Laplacian Ho is associated in the standard way 
with the quadratic form Q with domain Wq ' (Q). Similarly the Neumann Laplacian 
Hn is found from Q with the domain W^''^{Q). 

The usual Dirichlet and Neumann Laplacians (whose forms are defined on Wq'^{Q) 



N- 



and W^' (fi) respectively) will be denoted by —Ad and —A 
Formally, the differential operators are given by 

■^^ OXi \ OXi I 

i ^ ' 

and in particular for N = 2 

Hf = -a'^AEf 

although these are for illustration only and the rigorous definitions above are to be 
preferred. 

2.2. Manifolds with Fractal Boundary. The final section of [jl| gives the con- 
struction of a piecewise smooth homeomorphism between the ball equipped with 
a metric of the above type and S, the Koch snowfiake domain, equipped with the 
Euclidean metric such that the two metrics are Lipschitz equivalent (i.e. quasi- 
isometric). It is commented in this section that the technique seems to be restricted 
to two dimensions. This result allows one to establish a bound for the heat kernel 
K{t, X, y) of Q with the Riemannian metric, subject to Neumann conditions: 

< K{t, X, y) < kt-^ for < t < 1 

where k is some constant. This is possible because there is a body of knowledge 
about the snowfiake domain. More specifically it has the ly^'^(S') — > ly^'^(R^) 
extension property being a quasidisc (see e.g. Pj). 
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2.3. Other Examples. In [^ Pang considers singular elliptic operators in a sim- 
ilar situation. The operator acts in L'^{Q,ri^{x)d'^x) and is derived in the usual 
way from the quadratic form 

Qiif) = I r/lV,,/||d^a; where / G CT{n). 
Jn 

Q is again a smooth domain in M^, A and fi are constants and r/ is a smooth 
function such that ri{x) ~ cr{x) as x ^ dQ. The proofs are based on differential 
operator techniques and he systematically establishes a variety of L^ and L^ spec- 
tral properties such as conservation of probability and essential self-adjointness. 
The key point is that these properties depend only on the values of A and //. It 
should also be mentioned that as well as the above case he also considered the 
case Q = M^ and r^ ~ (1 + |a;p)2 and in this case he was later able to weaken the 
smoothness assumptions on rj in WM. For results related to this case see also 0. 



3. RiEMANNIAN DISTANCE TO THE BOUNDARY 

We denote the Riemannian distance between two points x,y & flhy d{x,y). The 
Riemannian distance can be calculated in two distinct ways which we recall below. 
The first allows us to compute upper bounds. 

Definition 3.1. A smooth curve w : [a, (3] ^ Q is said to be regular if w'{t) ^ 
for t G (a, /3). A continuous map w : [a, P] ^ Q is called piecewise regular if there 
exists a finite subdivision a = ao < ai < ■ ■ ■ < a^ = P such that w|[a,_i,ai] is regular 
for i = 1, . . . ,k. Piecewise regular curves will be known as admissible curves 

Definition 3.2. If the Riemannian arc length is given by 



ds := gij{x)dxidx 



3 



then an admissible curve w := [a,/?] — > Vt with coordinates w*(t) has length (i{w) 
defined by 



using the summation convention. Now we define d{x, y) by 

d{x,y) := mi{i{w) : w is an admissible curve and w{a) = x and w{P) = y}. 

It can be shown that 

d{x,y) = sup{|'?/'(x) — V'd/)! : V' is Lipschitz continuous and \Vip\ < 1} . 

This can give useful lower bounds. 

Also considered in [^ is the completion of Q and it gives formulae for the geodesies 
connecting two points Vi, V2 G dfl. These are curves in Q that intersect dQ at right 
angles. One further important piece of notation is the following: 

d{x) := dist(a;, 90) 
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for X G f2. As usual 

dist(x, (9^2) := inf d{x,y). 

An important result is the relationship between d and a. 

Theorem 3.3. Let Q and 7 be as in the basic model. For x & Q we have 



d{x) 



a^-^ix) 



1-7 
Proof Let ip be given by 

M^) = -, 

1-7 

A calculation then shows that in the sense of distributions 

|Vi,^o|| = ^"'lVi,a|| 
and 

iv^or = ^'^v^^oii < 1. 

Now, ipQ G W^'^ and so ipo is Lipschitz continuous, satisfies {Vipol < 1 and hence 

d{x) > 



a^-^/ix) 



1-7 

Next, given a point a; G f2, let A be the straight line segment that minimises the 
Euclidean distance to the boundary, i.e. the straight line segment that orthogonally 
intersects the boundary at the nearest point in the Euclidean sense. Choose the 
coordinate system so that the xi axis lies along A and the origin is at the point of 
intersection of A and dQ. In this coordinate system we have x = (a, 0, . . . , 0) for 
some a. 

Now for any point y = (t, 0, . . . , 0) lying along the straight line segment A we have 
(j{y) = t. Thus, the length of A is given by 

a^~^ a^-^ix) 



t-^dt 
1-7 1-7 



a^-^ix) 



Hence 

dix) < 

1-7 

and we have the result by combining the two inequalities. 
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4. The Hardy Inequality 

4.1. Generalised Hardy Inequalities. The Hardy Inequality is a useful tool 
when investigating the spectral properties of elliptic operators. It will play a role 
in what follows but first we stop to gather the essential definitions and results. 

We begin with M a Riemannian manifold, L a non-negative, second order elliptic 
operator acting in L^(M, dwo/) subject to Dirichlet boundary conditions and p a 
positive, continuous function defined on M which satisfies |Vp| < 1 at least in the 
distributional sense (p will be taken to be the distance to the boundary). Also let 
Q{f) be the quadratic form associated with L. 

Definition 4.1. We say that L satisfies a weak Hardy inequality with respect to p 
if there exists a constant c > and a constant a < oo such that 

|2 



[\IL^^ol<c'{Q{f)+a\\fr) 
Jo. P 



for all / in the domain of Q. The weak Hardy constant is then the infimum of all 
possible c such that this holds. 

If we may take a = in the above then we say L satisfies a strong Hardy Inequality, 
and the infimum is then called the strong Hardy constant. 

Further information may be found in the review by Davies ([0]) - also of note is the 
book by Opic and Kufner ([ly]) which establishes conditions for a Hardy Inequality 
to hold in a region with Holder class boundary. No attempt is made to quantify 
the constant in general. 

We conclude these introductory remarks with an important standard result. For a 
proof see . 

Theorem 4.2. IfVt C M.^ (usual Euclidean metric) is hounded with a C^ boundary 
and Ao is the associated Dirichlet Laplacian then there exists a &W such that 

in the sense of quadratic forms. 

4.1.1. HI for the Basic Model. We need the following result which is trivial to 
prove but an important observation. 

Lemma 4.2.1. 

if either side is finite. 

Proof Simply recall that dvol = a~^"'d^x and d{x) = "^ ^_ ■ 
Thus if we could establish a Hardy Inequality of the form 

(2) / M!^(2-iv)7d^^ < c' f / I V,,/||a(2-^)^d^x + a [ |/| V^^'^^^d^a; 
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then we would immediately have 

^\ol<c^(l-^?i I |VRf|la(2-^)^d^x + a' / I f I V'^^d^x 



n <P 



if { I I V^/||a(^-^)^d^x + a' / I/I V-^^d^ 



where a' = a||a^'^||oo and / G Wq' . 

To prove inequality (0) we work with the Hilbert space L^(r2, a^'^~-'^^"'d'^x) and use 
the following: 

Theorem 4.3. [^ Theorem 1.5.12] Suppose L is an elliptic operator on L'^{Q) 
and that there is a positive continuous function G Wj^^^ (Q) and a potential V in 
L^{Q) such that 

L(j)>V(l) 

in the distributional sense. Then we have 

L>V 

in the sense of quadratic forms. 

Now we can prove the first main theorem. 

Theorem 4.4. IfQ is convex (in the Euclidean sense) and equipped with the Rie- 
mannian metric then the strong Hardy Inequality holds with strong constant c sat- 
isfying 

c< ^(^-^) 
- 1 + {N-2)Y 

Proof Let L denote the Friedrichs extension of the operator given formally by 

Lf := -cr(^-2)^ y A f^i^-Nh^] 

^-^ OXi \ OXi I 

i ^ ' 

L is defined in terms of quadratic forms as usual. Now, let be defined by 

0:=a" 

where 

1 + (AT -2)7 

a := . 

2 

Then we have 

L0 = -a{a - 1 + (2 - iV)7)(T"-2|Vija|| - aa°-^Aija 
> -a{a - 1 + (2 - iV)7)a"~2 

in the sense of distributions. This calculation uses the fact that A^cr < for any 
convex set Vt. 
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Now, applying Lemma [4.3| we have 
and applying Lemma [4.2.1| gives 



Ja i(if " (l-(2-A')7Pin 

Theorem 4.5. T/ie wea/c Hardy Inequality holds for all Vt in the basic model with 
constant c satisfying 

2(1-7) 



c< 



1 + (AT -2)7' 



Proof Let $ be a smooth function that equals a near to the boundary. We now 
consider the differential operator K defined on C^ functions say, by 



Kf := -$(^~2h A (^$(2- 
dxi V 



dxi 



and again we denote by K its Friedrichs extension. 

The trial function is now chosen to be a smooth function which near the boundary 
satisfies 

= a" - a 

where a is as before. 

Then, 

Kcp _ -a{a - 1 + (2 - N)-f)a-^ - aa-^Asa + (2 - N)-fa'^-'' + a-'^Asa 

provided we are close enough to the boundary. 
Now observe that 



to conclude that 



Kl ^ (l + (iV- 2)7)^ ^,2 
(h - 4 



near the boundary. The behaviour of away from the boundary is not important 
for the existence of the weak Hardy Inequality and we conclude that 

Kct){x) > ('^ + (^-'^)^) a-\x)4>{x) - a0(x) 

for all X G fi provided a is chosen to be large enough. An application of Theorem 
01 establishes the final result. 
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Remark 4.6. The previous theorem indicates the general point that the weak 
Hardy Inequahty depends only on the local geometry of the boundary. To gain 
information about the strong version we must make assumptions about the global 
geometry. 

There is a link between the strong Hardy constant and the Minkowski dimension 
of the boundary. This link allows us to prove that the result given above is sharp. 

The following definition is taken from - it is not the usual definition but is 
sufficient for our purposes. 

Definition 4.7. We say that the boundary dQ has interior Minkowski dimension 
a if there exist positive constants Ci and C2 such that 

Ci£^-° < vol(5f],) < C2£^-" 

where 

dn,:={xen: dist(x, dn) < e} 

for small enough e > (vol(y4) denotes the Riemannian volume of the set A. Note 
also that the distance is calculated with respect to the Riemannian metric). 

Theorem 4.8. [|], Theorem 3.1] Let a denote the interior Minkowski dimension 
ofdQ. Then 

N ~1 <a = — ^— < N 

1-7 

The following result of Davies and Mandouvalos is now of interest: 

Theorem 4.9. [^ Theorem 6] If dQ has interior Minkowski dimension a, so that 
a > N — 1, then a, the strong Hardy constant of Q with respect to the Laplacian, 
satisfies 

c{2 + a-N)>2 

In conclusion: 

Theorem 4.10. Let Q be convex and let c denote the strong Hardy constant. Then 

,_ 2(1-7) 
1 + (AT -2)7' 



Proof Combine Theorems 4.4, 4.S and |4.9|. 



The range of values that c may take are of importance and we make a note of this 
now. 

Lemma 4.10.1. For 7 G [0, 1/N) the number 

c(7) = 2(l-7)/(l + (iV-2)7) 

satisfies 

1 < c<2 
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Proof Recall that 

c(7) = '^'-'^ 
^^' 1 + (iV - 2)7 

and < 7 < 1/N. 

A straightforward calculation shows that c'{'~f) < (where ' indicates differentiation 
with respect to 7), c(0) = 2 and c{l/N) = 2{N - 1)/{2N - 2) = 1. 

5. SoBOLEV Spaces and Norm Estimates 

There are two Hilbert Spaces which have featured thus far: the usual Euclidean 
space L^(n, d^x) which has norm 

ii/ii2,i.= /i/(^)rd^^ 

and the space L'^{VL,dvol) which is derived from the Riemannian metric and has 
the norm 



|/(x)pdt;o/= / |/(x)|V-^^d^x. 



Lemma 5.0.2. L^(i7,dfo/) C L^(i7,d^x) Moreover we have 

lk5^^||-i||f|UR< 



2,E 

Proof This follows from the fact that 

Turning to the quadratic forms we see that the usual Euclidean quadratic form 



QeU) = I |Vi./||d^x 
Jn 



and the Riemannian form 



QU) = I Wn'dvol = [ |V^/||a(^-^)-d^ 
Jn Jn 



X 



are equal when N = 2. Moreover, we can see from the previous lemma that 
QeU) < Q(/) in general. We now turn to the Sobolev spaces defined previously 
(see section ETll). We recall first a definition. 



Definition 5.1. We say fi is a regular domain if a generalised Hardy Inequality 
holds: i.e. there exist constants Ci,C2 with Ci > such that 

Hd > ci/(f - C2 

where Hd is the natural Dirichlet Laplacian and d the distance to the boundary. 



Remark 5.2. Theorems [4.2| and [4.5| imply that fi is a regular domain with both 



the Euclidean and the Riemannian metric. 



12 C. MASON 

The following theorem gives a complete description of the Sobolev space Wq ' . 
Theorem 5.3. If Q is a regular domain then 



w^'\n) = w'^\n)n^f:JJ-^ 



< oo 



Proof This theorem is essentially proved in |^, Theorem 1.5.6]. The proof involves 
a regularised distance function. 

Lemma 5.3.1. If f e Wq'e{^) then f e L'^{Q,dvol). Moreover, 

= \w-'-''yhE<\\fcr-%,E\w'-'-'''\ 



\oo- 



Proof If / G W^o'eI^) then Theorem pT3| ensures that 

-d^x < oo 



l/P 



a 



2 



i.e. a ^fE L^(f2, d'^x). Also a^ 2^t g -L°° and so Holder's Inequality immediately 
gives us that 

^-'"'^'fhE<\\f^-%,E\W'-^'''\ 



\cr 



oo- 



Lemma 5.3.2. For arbitrary dimension N we have 

and hence W^'^n) C W^^U^). 
Moreover, if N = 2 we have 

Proof The first two statements follow immediately from the discussion about the 
norms and the quadratic forms. Now, if / € Wq'^{Q) then / G L'^{Q,dvol) by 
previous lemma. Next observe Q is independent of 7 when N = 2 and so it follows 
that / G iy^'^(r2). Finally we observe that 



n " Jn (^ 



-dvol = (1 — 7) / — — d X < 00 



and thus we conclude that / G Wq' (i7). 

6. Boundary Perturbation 

We now return to the boundary perturbation problem as described in the introduc- 
tion. We recall that in the case that a weak Hardy Inequality holds in the domain 
with Hardy constant c > 2 the result is as follows: 

The difficulty in applying these results to our situation is that 1 < c < 2 and 
indeed c < 2 whenever 7 > 0. Of course, the Hardy Inequality remains true with 
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the larger constant c = 2 and we could simply apply these results to get the fact 
that 

An,e- A„ = 0(e). 

This is less than satisfactory. Numerical experiments and specific examples suggest 
a rate of convergence which is faster than this - in the following section we set up 
an example that may be used to show a faster rate. 

The results in this section will be general in character although they will ultimately 
be applied to the basic model. For this reason we think of Q as being some 
Riemannian manifold with Dirichlet Laplacian H^ and Hardy Inequality 

where d is the distance to the boundary and c > 1 some constant. 

We will modify the approach in [^ . This uses the original eigenvectors {/)j of H^ 
to approximate the eigenvectors of H^jy in the variational formula. We define a 
rapidly decreasing function /i : fi — *• [0, cxo) by 

r if < d{x) < e 

^{x):= } ^fc iie<dix)<2e 
[^ 1 otherwise. 

We note that < /i < 1, |V;u| < - and /i has support in Q^- 

The effect of applying this cut-off function is captured in the following theorem. 

Theorem 6.1. For f EDomfHi)) the following holds 

1. fif eDom(Q), 

2. QM < Qif) + 2 /^..,<,(,.)<,,^ \Vf\Mvol + 2e~^ 4^,<,(,.)<,,^ Ifl'dvol, 



Proof Theorem |6.1| : 1 More generally we have the following : 

Lemma 6.1.1. P, Lemma 2] // / EDom(Q) and w G W^'^ then wf EDom(Q) 

Proof [Lemma |6.L1 1 



Q{wf)= / {\wWf + fWwY)dvol 
Jn 

< [ {2w^\Vf\^ + 2\ff\Vw\^)dvol 



<2\\w\\lQ{f)+2\\Vw\ 



Now we may give the conclusion of Theorem |6.1[ 1. 

It has already been noted that fj, G W^'°° and so the lemma can be applied. 
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Proof Theorem |6.1| :2 Let S := {x : e < d(x) < 2e}. Then we have 



Q(/i/)-W)< / \vifif)\'dvoi 

's 



<2 / /x'|V/|Mwo/ + 2 / \Vfi\'\f\'dvol 
's Js 



<2 



/iv/r 

Js 



dvol + 2e-^ / \f\Mvol 



Proof Theorem 16.11:3 



ll/^/l|2p<||/-/^/||; 



'l-^f\f\Mvol 



n 



< 



'dvol. 



' {x:0<d{x)<2e} 

This completes the proof. 

The key point now is to obtain boundary decay estimates by considering the fol- 
lowing integral 



J{x 



\ffdvol 



' {x:d{x)<£} 

making only the assumption that a Hardy Inequality of the form 

d-^ <c\HD + a) 

holds in the sense of quadratic forms for c > 1. 

Although we have relaxed the condition on c this is not without incurring some 
further cost. The estimates achieved in |^ were valid for functions lying in the 
domain of the operator. The following theorem captures the key results of ||^. 

Theorem 6.2. Let f eDom(HD), c^{Hd + a) > rf-^ j^r c > 2. Then 



L 



\f\'dvol<^+'/^\\{HD + a)fU{HD + af'^fh 



{x:d{x)<e} 



and 



L 



\Vf\^dvol<{c^''^ + ^'%l + cf+^/^)\\{Ho + a)fU{HD 



vl/c 



fh 



{x:d(x)<e} 



Proof See |^, Theorems 4 and 8] 

Our results will be even more restrictive in that they demand that the functions lie 
in the space Dom((ilf£)+a)^/^+^/^). This is not a serious restriction when the results 
are applied to the perturbation of domain problem. In this we are only concerned 
with the eigenvectors of H£, which lie in Dom((if£) + a)*) for all s > 0. Indeed, 
when the results are extended to more general elliptic operators with measurable 
coefficients the distinction becomes perhaps even less important because in this 



PERTURBATION OF DOMAIN: SINGULAR RIEMANNIAN METRICS 15 

case it is extremely difficult to determine even Doin{Ho) and certainly it cannot 
usually be identified with one of the standard Sobolev spaces. 

The following result is taken from , it does not depend critically on the value of 
c and therefore is also applicable in our situation. 

Definition 6.3. We define the function w by 

w{x) := {ma.x{d{x),e})~^^'^. 

Lemma 6.3.1. [^ Lemma 3] // / (zDom(H£,) then 

I ( J^ - l^^l') l/l'd^«^ ^ {Hf,w'f) + a\\wf\\l 



Proof By the Hardy Inequality and Lemma |6.1.1| we have wf GDom(Q) and 

Jn 



— ^dt;o/ < Q{wf) + a||M;/||2. 



Moreover, 



Qiwf) = {Hnf,w'f) + [ \Vw\'\f\'dvoL 

Jn 

This calculation may be found in slightly more detail in but is routine. 

The main theorem in this section differs from a similar result in P] in that we can 
ignore the restriction that c > 2. 

Theorem 6.4. Let f G Dom{{HD + a)^/^+^/') where 1 < c < 2. Then 

\{Hnf,w'f)+a\\wf\\l\<c'/^UH + aY/'+'/^fh\\iHn + aY/'fh 

Proof We have 

in the sense of quadratic forms. Thus 

{Hd + ay^w^^HD + a)-^2 < c^ 
and so 



1 



Ww^iHo + a) 2 II < c. 

We also have 

(3) w^^^-'^ < (d-2)^ < {c'iHD + a))"-^ 

since < (2 — c)/c < 1. This uses the fact that for non-negative self adjoint 
operators A and B the quadratic form inequality < A < B implies < A" < 5" 
for all < a < 1. It is at this point that the assumption c > 1 is important. For 
a proof of this fact see 0, Lemma 4.20]. 

Now from (B) we also have the fact that 



2-c ,, 2-c 



w'-'^iHD + ay—W <c 
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Finally we make the definition Ha '■= H^, + a and then 

\{Hnf,w'f) + a\\wf\\l\ = \{{H + a)f,w'f)\ 

= \{w''~'h;^ H^ hj,w'h;~'hI f)\ 

The theorem can now be applied in the same manner as . The following theorem 
is Theorem 4 from |^ using the above estimate in the proof. 

Theorem 6.5. Let f E Dom{{HD + 0)^/^+^/"). Then 



{x:d{x)<£} 

Proof Combining Theorems |6.3.1| and |6.4| we have 



in \c'(P 
If d{x) > e then 



Thus 



Also, if d{x) < e then 



Hence 



\Vw? < c-'d-'~'/' ""' 



c^d^ 



w 



2 



- ivwr > 



c^d^ 



C^d^ ' ' C%2/C(i2- 



J{x:d(x)<e} " Jn \C^d J 

The result follows immediately. 

The previous results allow us to obtain similar estimates for V/. 

Theorem 6.6. Let f G Dom{{HD + ay/^+^/^). Then 



j 

J\x 



{x:d{x)<£} 

where Ci := c^l^ + c^'^{l + cf+'^'^. 
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Proof This is similar to Theorem 8 in 0] . Again we use Theorem |6]J in place of 
Theorem 1 of 0] at the appropriate point. 

As previously mentioned, when applying the results achieved in the last section 
we consider only the eigenf unctions of the original operator Ho and use them to 
approximate the eigenfunctions of the restricted operator Hcfr. 

Definition 6.7. Let 0i, 02 • • • 0n denote the eigenfunctions oi Hd, where {(pi, (pj) = 
6ij, corresponding to the eigenvalues Ai, A2, . . . , A^. 

Now, define the n-dimensional subspace L^ by 

Ln := span{0i,...,0„} 



We now state a corollary to Theorems |6.5| and |6.6| 
Corollary 6.8. Let f e L^. Then 



I \j I Liuui, ^ o c y/\n -r "; \\J 112 

J {x:d{x)<£} 



\f\'dV0l < c2+2/%2+2/c^^^ _^ ^)l+l/c|| f ||2 



and 



I \Vf\^dvol<Cie^''{Xn + a) 

Jfx:d(x)<£\ 



' {x:d{x)<£} 



Proof It is immediate from the spectral theorem that 

||(i7z. + ar/||2<(A„ + ani/||2 
for f E Ln. Now apply Theorems |6]^ and |0| . 
Theorem 6.9. Let c > 1. Then, there exist constants Cn such that 

^n,£ ^n — "^n^ 

for £ > small enough. 

Proof We use the variational estimate 

K,e < sup ^ . 

Combining Theorem |6.1| and Corollary ^]8| we find that 

^ ^ A„ + c-£^/-(A„ + a)i+V- 
"'' " 1 - ci+i/-ei+i/-(A„ + a)5(i+Vc) 
where c' = 2(ci + c^+^Z'^). 

Corollary 6.10. Let Q be as given in the basic model, satisfying a Hardy Inequality 
with constant 

,_ 2(1-7) 
l + (A^-2)7' 

Then there exist constants Cn such that 

\ _ A < r ^(i+{A'-2)7)/(i-7) 
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for e > small enough. In particular if N = 2 we have the following 

K,e - An = 0{£^^) as £ -^ + . 

7. The Rotationally Invariant Case 

We now consider a specific example which is particularly tractable and may be 
used to show that the power achieved in the previous general result is sharp. We 
use the unit disc as the underlying region, i.e. 

Q:= {x: \x\ < 1} C M^ 

and so the Riemannian metric is given by 

ds^ = (1- \x\)-^^\dx\^ 0<7< i . 

Using the orthogonal group we make the decomposition 

where 

L V 2n 

Each subspace L^ is invariant with respect to the Laplacian and so we can restrict 
our attention to each subspace. This uses the well known process of separation of 
variables. Thus we can consider only 

"1 
2_ / \g(^r)\'^(l — r)~'^"'rdr 



Qnif) 



/ (Wir)\^ + ^\9{r)nrdr. 



We reiterate that the case N = 2 causes an important simplification: the quadratic 
form and its domain are independent of 7. Formally, the associated differential 
operator is given by 

(4) Lg := -(1 - rf (i|- (r|) - ^g^ on ^^((0, 1), r(l - r)-^-rfr). 

Next we reformulate the problem is terms of a Schrodinger operator. 

Lemma 7.0.1. Under suitable transformations we have the following expressions 
for the norm and quadratic form: 



l^= / \Ht)\'dt 

Jo 

Q{h)= r [\h\t)\^+v{t)h{t))dt 
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where 

1 
a 



1-7 
and 



4(1 -jf {a -tr V ij (l-il-i)")''' 

Proof Let x = w{t) = 1 — ((1 — 7) (a — t))" so that t = corresponds to a; = and 
t = a corresponds to a; = 1. We work with g G C^{Q) and then prove the final 
result by approximation. The norm and form become 



where 



\\9\\l= / \9{w{t))\a{t)dt 
Jo 

QU)= r {\g\w{t))\^dt+\g{w{t))\b{t))a{t)dt 



wit) 
a{t) :-- ^ ' 



w'{t) 



b(t) := n 



w{t) 



2 



Make the substitution h{t) := g{w{t))^Ja{t), then 



h = ga2 +-g— 

,1 la' 
= ga'2 +-h—. 
2 a 



Thus 



r\gfadt= r (\h'?+-i-f\hAdt- rhh'-dt 

io V A a J Jq a 



/|2 I -"- /" \2|L|2 \ ^j. I / -"-1^12 






Immediately we have 
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A routine, but messy, calculation establishes that in terms of w we have 

^ ' A\w' ) 2w' \ ^) \w 

The final form of the potential is found after utilising the formulae: 



t^" 



w{t) =1 - ( 1 - 

urn ^ (i - ^) 



a 



«;"(t) = - 7 (l - ^y ' 

w"'{t) =7(27 - 1) (l - ^' " 

After some calculations we see that we can reduce the problem to considering 

(5) -h"{t) + V{t)h{t) = Xh{t) 

for 

< t < a 
and 

<t < a- e 
where 

v(t) = "" ~y^ + Oir^) as t -^ + 

^<" = Z'-tyHc-t)^ + ''«" - oft) as t ^ a - . 

In another paper, 0, we deal with general differential equations of this type. The 
techniques and proofs are different from anything in this paper. We refer to |^ for 
details and here simply quote the final result as it applies to our model. 

Theorem 7.1. Let X„ be the nth eigenvalue of equation (j^ on the full interval 
(0, a) and A„_£ the nth eigenvalue on (0, a — e). Then 

A„,e = A„ + Cn€^=^ + 0(£i^) OS £ ^ 0+ 

for some constant Cn- 



The power in this expansion agrees with that in Corollary |6.10 and shows that it 
gives the best result possible. 
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